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Heavy nuclei bombarded with protons and deuterons in the 1 GeV range have a large probability 
of undergoing a process of evaporation and fission; less frequently, the prompt emission of few 
intermediate-mass fragments can also be observed. 

We employ a recently developed microscopic approach, based on the Boltzmann-Langevin trans¬ 
port equation, to investigate the role of mean-field dynamics and phase-space fluctuations in these 
reactions. 

We find that the formation of few IMF’s can be confused with asymmetric fission when relying on 
yield observables, but it can not be assimilated to the statistical decay of a compound nucleus when 
analysing the dynamics and kinematic observables: it can be described as a fragmentation process 
initiated by phase-space fluctuations, and successively frustrated by the mean-field resilience. As an 
extreme situation, which corresponds to non-negligible probability, the number of fragments in the 
exit channel reduces to two, so that fission-like events are obtained by re-aggregation processes. 

This interpretation, inspired by nuclear-spallation experiments, can be generalised to heavy-ion 
collisions from Fermi to relativistic energies, for situations when the system is closely approaching 
the fragmentation threshold. 


I. CONTEXT 

Several decades passed from Serber’s early descrip¬ 
tion [1] of nuclear reactions induced by nucleons and light 
nuclei at few hundred MeV per nucleon. In their stan¬ 
dard outline, such reactions, generally called spallation, 
could be described as a fast excitation of an atomic nu¬ 
cleus, followed by a slower decay process W- depending 
on the phase space available HIS], the system undergoes 
a sequence of more or less asymmetric splits 00 rang¬ 
ing from particle evaporation to fission. It was found 
already in pioneering studies 0E], that the most ex¬ 
cited systems can also produce intermediate-mass frag¬ 
ments (IMF), and lead to a richer phenomenology com¬ 
parable with nucleus-nucleus collisions and nuclear frag¬ 
mentation [IHHII]. Further research focused on the study 
of thermodynamic observables from spallation reactions 
in the relativistic domain |22H25j , in connection with the 
liquid-gas phase transition in nuclear matter [26l[27], and 
in parallel with the research on the multifragmentation 
process observed in ion-ion collisions in the Fermi-energy 
domain [281482] . 

Several fields of application, from energy and environ¬ 
ment to neutron sources and exotic beams, stimulated 
intense research on protons and deuterons in the 1 GeV 
range impinging on heavy nuclei. The production of some 
specific light nuclides with large kinetic energy resulted of 
great relevance in several technical issues (radiation dam¬ 
ages, fragilisation of structural materials in accelerator- 
driven systems, side effects in medical hadron-therapies). 
This despite the minor contribution of the whole IMF 
production to the total reaction cross section, which was 
found to amount to few millibarn. In more recent ex¬ 
periments, the possibility of correlating isotopic cross 
sections to high-resolution kinematic observables allowed 
tracking the process of the IMF production: it was ad¬ 


vanced that it could originate from a melange of two pro¬ 
cesses, multifragmentation from the most excited config¬ 
urations |83j . and asymmetric fission [34] . 

This work proposes a fully dynamical description of 
the process, within a Boltzmann-Langevin transport ap¬ 
proach, with the aim of probing the mechanism which 
rules the physical process, through the analysis of kine¬ 
matic and correlation observables. 


II. THE QUEST ON THE ORIGIN OF IMF 
IN SPALLATION 

The IMF production in spallation, especially in the lA 
GeV range, is a process at the threshold between multi¬ 
fragmentation and compound-nucleus decay. To charac¬ 
terise the process, an ideal experimental approach should 
measure event-by-event particle-particle and kinematic 
correlation observables at high-resolution. This goal has 
been only partially achieved, so that the lack of resolu¬ 
tion or the missing of some correlations opens the way 
for a variety of physical interpretations, ranging from at¬ 
tributing all IMF’s to sequential fission processes to the 
opposite extreme that all IMF’s signal multifragmenta¬ 
tion events. In the following we will refer to two specific 
experimental approaches, implying that the same consid¬ 
erations could follow from other experimental analyses. 


A. Experimental observables: an example 

The first experimental approach we focus on is a cam¬ 
paign of inverse-kinematics experiments, performed at 
GSI (Darmstadt): the FRagment Separator [35] (FRS), a 
high-resolution magnetic achromat [36] , was employed to 
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measure the nuclide production in several spallation re¬ 
actions at around 1 AGeV [33l IHTMO] . In some systems it 
was possible to measure the isotopic cross sections of the 
IMF and the corresponding zero-angle invariant-velocity 
distributions; these distributions are constructed by se¬ 
lecting only the velocity vectors aligned along the beam 
direction (they are evidently different from longitudinal 
projections of the whole velocity distributions). Fig. 
presents some of those experimental results taken from 
ref. [49] , and resumes few essential steps of the data anal¬ 
ysis [5T] . 

As evident in Fig. in the IMF region the veloc¬ 
ity spectra of a given nuclide is the composition of two 
shapes: one is concave (showing two peaks), the other 
is convex. The two contributions were also found to 
have different relative shares as a function of the nuclide. 
Though, the presence of two different contributions is ev¬ 
ident only in the kinematics, while the nuclide production 
selected for either the concave or the convex mode con¬ 
tribute to the same region of the nuclide chart, in general 
situated in the neutron-rich side with respect to /3 sta¬ 
bility: this is the reason why the presence of these two 
modes was invisible in many experiments. Moreover, the 



FIG. 1. (a) Isotopic cross sections a measured in the IMF re¬ 
gion for the system '^^^Xe-l-p at lAGeV, from the experiment 
of ref. |49 |. (b) Zero-angle invariant velocity distributions 

CTi(uf|)/CT extracted from the experimental data for some iso¬ 
topes [^. (c) Analysis detail for one isotope, ^°F: invariant 
velocity plot on a plane containing the longitudinal axis. Two 
contributions are present with a recoil mismatch: a Goulomb 
ring with radius Upeak and a gaussian-like distribution, (d) 
Reduction of the two-dimensional plot to a distribution along 
the beam axis: two kinematic contributions to the invariant- 
velocity spectrum appear as a convex mode (Coulomb ring) 
and a concave mode (gaussian-like), respectively. 


two contributions are associated with a shift in the mean 
value of the spectra, indicating incompatible values for 
the mean momentum transfer: this latter reveals in fact 
the violence of the entrance channel. 

The concave shape is reflected in a Coulomb-shell ve¬ 
locity distribution probed at zero angle. The radius of 
the shell is compatible with a fission kinematics [^ and 
the mean value evolves coherently with empirical system- 
atics for the mean recoil momentum as a function of the 
fissioning system (i.e. compatible with the systematics 
of Morrissay [55]). 

The convex shape is one broad hump, often asym¬ 
metric, which signs the folding of many different contri¬ 
butions; it is associated with a mean recoil momentum 
which does not follow any empirical systematics. This in¬ 
dicates that the fissioning configuration is not achieved, 
either because the multiplicity of fragments is not equal 
to two, or because the kinematics is not consistent with 
a conventional fission configuration. In this respect, the 
convex shape is rather compatible with multifragmenta¬ 
tion. 

Such observations led to the conclusion that the IMF 
production should combine asymmetric-fission (concave 
shapes) and multifragmentation (convex shapes). On 
the one hand, this description was rigorously established 
because these two velocity contributions were measured 
at the same time for each single nuclide m- On the 
other hand, the problem of assuring such interpretation 
was that particle-particle correlations and IMF multiplic¬ 
ities were not measured in the inclusive approach. Espe¬ 
cially for the convex shape this information is important 
as multifragmentation is usually associated with a large 
number of IMF of similar size. 

A second experimental technique was adapted to ob¬ 
tain this information. Some of the systems previously 
measured inclusively at the FRagment Separator were 
successively measured again with an exclusive approach 
in the Spalladin experimental campaign |51H55] at GSI 
(Darmstadt). These experiments indicated that IMF are 
observed in events where fragment multiplicity is preva¬ 
lently equal to two, and that events with larger multiplic¬ 
ity were more rare. On the one hand, this confirmed that 
there are actually two contributions to the IMF produc¬ 
tion, a binary channel and a higher multiplicity channel. 
On the other hand, the events exceeding two IMF were 
not easily identifiable with ordinary multifragmentation 
because of the low IMF multiplicity and the size asym¬ 
metry. This encouraged interpretations fully relying on 
statistical models, where the IMF production is obtained 
either from a sequence of (asymmetric) fission contribu¬ 
tions, or sampled from an ensemble of possible multifrag- 
mented configurations. In general, these approaches are 
both an efficient workaround because they define directly 
the outcome of the reaction on the basis of the involved 
excitation energy, and they can yield quite similar results 
for the IMF production despite implying different phys¬ 
ical pictures. Then the question arises whether we can 
include spallation in the multifragmentation picture. 
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It is our intention to illustrate that the key to achieve 
a fully coherent understanding of this apparently self¬ 
contradictory experimental information, conciliating the 
possibility to emit only few IMF’s with the observation 
of new kinematic properties, is a microscopic dynamical 
description of the process. 


B. Need of a dynamical description 
to address the quest 

Avoiding prominent mechanical contributions from the 
entrance channel which characterise ion-ion collisions at 
Fermi energies, spallation is a favoured processes to pro¬ 
duce a thermalised remnant. It is therefore well estab¬ 
lished that the usual hypothesis of a hot and fully equi¬ 
librated source is perfectly adapted to spallation and in 
fact statistical models proceeding from this assumption 
are well suited to reproduce large part of the experimen¬ 
tal observables [33]. On the other hand, it is not the 
purpose of these approaches to investigate the reaction 
mechanism which leads to the final configuration. 

To directly address the quest, we proceed through a 
dynamical description avoiding any a priori assumption 
on the degree of equilibration of the system at a given 
reaction time. 

In particular, the dynamical approach allows to de¬ 
scribe the possible onset of mechanical instabilities in a 
fermionic system: this is a general process which char¬ 
acterises Fermi liquids at low densities and which results 
in inhomogeneous density patterns. As we will explore 
in the following, the incident energy of a light projectile 
can actually be large enough to produce the thermody¬ 
namical and density conditions for unstable modes to get 
amplified, i.e. to enter the nuclear spinodal region |57j . 
The dispersion relation [58j defines how the growth rate 
of these unstable modes connects to the mean-field po¬ 
tential. In this direction, a one-body approach based 
on an efficient description of the dispersion relation is 
well adapted, and a stochastic treatment is necessary 
to sample the variety of possible dynamical trajectories 
that unstable conditions may produce. With the purpose 
of investigating fragment formation, stochastic one-body 
approaches have already been applied to spallation [53] . 
but with incomplete success because fluctuations were 
not treated in full phase-space. As an improvement, in 
this work we employ the Boltzmann-Langevin One Body 
(BLOB) approach jSD] which has been constructed un¬ 
der the explicit constraint of describing the fluctuation 
mechanism in full phase-space by solving the Boltzmann- 
Langevin (BL) transport equation in three dimensions. 
The model was firstly applied to dissipative central ion- 
ion collisions [sniinii at Fermi energies, to investigate the 
transition from incomplete fusion to multifragmentation 
in proximity of the low-energy threshold for multifrag¬ 
mentation, which we expect spallation can also probe in 
the lAGeV regime. 


III. BOLTZMANN-LANGEVIN ONE BODY 

DESCRIPTION OF A SPALLATION SYSTEM 

Within the BLOB transport model we follow the dy¬ 
namics of the target nucleus, after the interaction with 
the light projectile. We describe the V-body system with 
a one-body Hamiltonian H supplemented by a fluctu¬ 
ating contribution to account for the unknown V-body 
correlations. The BL equation describes the time evolu¬ 
tion of the semiclassical one-body distribution function 
/(r,p,t) in its own self-consistent mean field: 

dtf-{H[flf} = i[f]+5I[f]. (1) 

The left-hand side gives the Vlasov evolution for / and 
the right-hand side introduces the residual interaction, 
which also carries the unknown N-body correlations. 
This latter contains the average Boltzmann hard two- 
body collision integral /[/] and the fluctuating term of 
Markovian type SI[f], also written in terms of the one- 
body distribution function [62] . 

The propagation of the one-body distribution function 
is described through the test particle method and em¬ 
ploys a Skyrme-like (SKM*) effective interaction [H3] . 
defined according to a soft isoscalar equation of state (of 
compressibility K = 200 MeV); the potential energy per 
nucleon E-pot/A is defined as 

( 2 ) 

with u = p/po and /3 = {pn — Pp)/p, where po and 
{pn — Pp) are the saturation and the isovector density, re¬ 
spectively. The parameters are A = —356 MeV, 5 = 303 
MeV, CT= 7/6. Surface effects are accounted for by consid¬ 
ering finite width wave packets for the test particles em¬ 
ployed in the numerical resolution of Eq.(l). The explicit 
term added to the potential energy is tuned to reproduce 
the surface energy of nuclei in the ground state [33]. A 
linear (stiff) density dependence of Epot is considered by 
choosing Cgymip) = constant = 32MeV. 

The fluctuating term SI[f] acts as a dissipating force 
during the whole temporal evolution of the process and 
introduces fluctuations by exploiting V-body correla¬ 
tions. In the BLOB procedure, binary collisions involve 
extended phase-space agglomerates of test particles of 
equal isospin A= oi, 02 ,..., B= 61 , 62 ,... to simulate nu¬ 
cleon wave packets: 

/[/] + SI[f] =9 /clH IF(ab«cd) F(ab^cd) , 

(3) 

where g is the degeneracy factor, W is the transition rate, 
in terms of relative velocity between the two colliding ag¬ 
glomerates and differential nucleon-nucleon cross section 

IV(AB.r^CD) = I ?;a - UBI ^ , (4) 

and E contains the products of occupancies and vacancies 
of initial and final states calculated for the test-particle 
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agglomerates 


_F(ab->cd) 


( 1-/a ) ( 1-/b )/c/d-/a/b( 1-/c ) ( 1-/d ) 


Since Ntest test particles are involved in one collision, 
and since those test particles could be sorted again in 
new agglomerates to attempt new collisions in the same 
interval of time as far as the collision is not successful, the 
nucleon-nucleon cross section contained in the transition 
rate W should be divided by Ntest- cf = cTNN/fVtest- In 
this work the unn is taken equal to the free nucleon- 
nucleon cross section, with a cutoff at 100 mb. Moreover, 
the differential cross section depends on the scattering 
angle according to the prescription of ref. [M] . 

All test-particle agglomerates are redefined at succes¬ 
sive intervals of time in phase-space cells of volume h?] 
in their initial state they correspond to the most com¬ 
pact configuration in the phase-space metrics which does 
neither violate Pauli blocking in the initial and in the 
final states, nor energy conservation in the scattering. 
The metrics of the test particle agglomerates is defined 
in such a way that the packet width in coordinate space is 
the closest to where corresponds 

to the screened cross section prescription proposed by 
Danielewicz [65j . which was found to describe recent ex¬ 
perimental data [66) . In this way, the spatial extension 
of the packet decreases as the nucleon density increases. 
The nucleon-nucleon correlations produced through this 
approach are then exploited within a stochastic proce¬ 
dure, which consists in confronting the effective collision 
probability W y- F with a random number. When the 
scattering is successful, a precise shape-modulation tech¬ 
nique m is applied to ensure that the occupancy distri¬ 
bution does not exceed unity in any phase-space location 
in the final states. Such constraint avoids that the Pauli 
blocking could be violated, and it imposes to pay special 
attention to the metrics of the phase space (see discussion 
in ref. [HH]!. 

As a consequence, fluctuations develop spontaneously 
in large portions of phase space, with an amplitude vari¬ 
ance equal to /(I — /), at equilibrium, in a phase-space 
cell of volume h^. This leads to a correct Fermi statistics 
for the distribution function /, in terms of mean value 
and variance. Calculations in a periodic box for unstable 
nuclear matter, in one dimension [69] and in three dimen¬ 
sions m have shown that the BLOB approach describes 
the growth rate of the corresponding (spinodal) unstable 
modes, related to the form of the mean-field potential, as 
ruled by the dispersion relation. Thus the BLOB model 
for heavy-ion collisions is constructed as based on the 
efficient description of these aspects. 


tion. As usual, the system is initially defined by organis¬ 
ing the test particles in a minimum-energy configuration 
in accordance with the form chosen for the nuclear inter¬ 
action. In order to define the heated system, this config¬ 
uration is redefined by processing a simplified cascade in¬ 
duced by the incoming light projectile: the amount of en¬ 
ergy deposited by the projectile in traversing the nucleus 
is calculated as well as the corresponding distribution in 
phase space. A time-dependent calculation would require 
very small time steps and a relativistic formalisation of 
the dynamics. Due to the rapidity of the spallation pro¬ 
cess with respect to the dynamics of the heated system, 
it is convenient to reduce the cascade to an approximate 
description where only test-particles from the incoming 
relativistic projectile are followed along space trajectories 
and target test-particles are not displaced during the cas¬ 
cade. In practice, this simplification is made by reducing 
the cascade process to a calculation of the energy loss 
of the projectile, modifying the momentum landscape of 
the target test-particles without processing any time evo¬ 
lution of the system in coordinate space. For relativistic 
projectiles this choice is not incompatible with the obser¬ 
vation that the projectile leaves the target nucleus before 
that the swarm of the first fastest ejectiles appears at the 
surface of the target nucleus [I11IZ2]- 

Fig. [2] shows a bunch of spacial cascade paths cor¬ 
responding to a ^°®Pb target nucleus bombarded by 1 
GeV proton projectiles with a central impact parameter; 
the resulting excitation energy distribution corresponds 
therefore to the most violent events. The test particles 
composing the projectile hit the nucleus along the pro¬ 
jectile direction within a cylinder of radius equal to the 
projectile radius. Each projectile test particle produces 
a cascade path inside of the target nucleus, which is re¬ 
defined at each scattering occurrence: after scattering, 
the projectile test-particle trajectory continues to be fol¬ 
lowed along the fastest scattered particle and the other 



A. Definition of the heated system 

Differently from heavy-ion collisions at Fermi energy, 
in this application to spallation only the dynamics of the 
heated heavy nucleus is followed after a suitable initialisa- 


FIG. 2. Energy-deposition map with one bunch of spacial 
cascade trajectories in relief, calculated for a central impact 
parameters in the reaction p-|-^°®Pb at lAGeV. Each trajec¬ 
tory is associated to one test particle of the incoming projec¬ 
tile. 
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particle, after being assigned a new momentum, is no 
more followed. Between two scattering points the path 
is a straight segment. All cascade paths traced by the 
projectile test particles are followed in coordinate space 
simultaneously. For a couple of target and projectile 
test particles, the collision is searched according to the 
closest-approach criterion applied to the corresponding 
centre-of-mass energy y/s |64j and by using the same 
nucleon-nucleon cross section used for the transport cal¬ 
culation. All collisions are considered as elastic scatter¬ 
ings; the model could be improved by including the A 
production-absorption mechanism, but we consider the 
present simplified treatment sufficient for the purpose of 
obtaining the excitation energy of the target nucleus. A 
strict Pauli-blocking condition here is imposed, so that 
only scattering events which create a hole and a particle 
outside of the Fermi sphere are accepted; otherwise, the 
target test particle could participate to a scattering with 
another projectile test-particle. When the cascade tra¬ 
jectories hit the inner potential boundary of the system, 
they can traverse the boundary according to the corre¬ 
sponding transmission probability, calculated with the 
relativistic formalisation proposed in ref. |73j ; the poten¬ 
tial depth used for calculating the transmission probabil¬ 
ity is 40MeV, which represents the average value charac¬ 
terising the bulk of the system. This transmission prob¬ 
ability is used to calculate an additional portion of the 
total energy of the projectile, which is considered dissi¬ 
pated in the target system and which corresponds to the 
reflected wave. 

While the coordinate space is frozen to its initial con¬ 
figuration, the initial momentum distribution is updated 
according to the cascade scatterings. The energy de¬ 
posited by the projectile in the system is then obtained 
by considering the momentum variation, supplemented 
by the reflection contribution at the potential boundary 
of the system. Accordingly, as shown in fig. for the 
reaction p-|-^°®Pb at 1 GeV, the initial integrated kinetic 
energy distribution Pq{E\^) is modified into a new distri¬ 
bution Pg(Ak). 



FIG. 3. Effect of the cascade in modifying the momentum 
space from the initial state configuration (dashed line) to the 
excited configuration (full line), calculated for one event in 
the reaction p-|-^°®Pb at 1 AGeV (see text). The panel on 
the left represents the nucleon energy distribution, whereas 
on the right the integrated energy distribution is presented. 




FIG. 4. Upper panel. Evolution of the mean number of emit¬ 
ted nucleons per interval of time in the reactions p-|-^''®Pb at 
1 GeV. Lower panel. Evolution of the mean excitation energy 
calculated for the fraction of bound matter in the reaction 
p_l_ 208 p^ Central, intermediate and peripheral impact pa¬ 
rameters are tested. The width of the bands give the standard 
deviation around the trajectories for the central collisions; 
other trajectories have a comparable standard deviation (not 
indicated). 


The cooling process of the excited system is then fol¬ 
lowed in time with BLOB. The latter is shown in hg. 
setting the reduced impact parameter 5red (impact pa¬ 
rameter divided by the target radius) equal to 0, 0.5 
and 1 for the system p-|-^°®Pb at 1 GeV. The evolu¬ 
tion of the mean fraction of bound matter (dAbound)/dt 
tracks the mean number of emitted nucleons per inter¬ 
val of time: central and intermediate impact parameters 
act almost equally in removing a large part of nucleons, 
while peripheral collisions favour the formation of heavier 
remnants. The corresponding information is carried by 
the evolution of the mean excitation energy per nucleon 
{E*/A) averaged over all portions of bound matter in the 
system. 


IV. NUCLIDE PRODUCTION 
AND KINEMATICS 

The model described above was applied to six systems, 
chosen because close to some signihcant experimental 
data and because they constitute a series of successive 
variations of only one parameter among projectile, tar- 
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get and energy: ^°®Pb+p at 1 A GeV, at 750 

A MeV, 208pb+d at 750 A MeV, at 750 A 

MeV, ^^®Xe+p at 1 A GeV and ^^^Xe+p at 1 A GeV. 
The dynamical calculations were performed reducing to 
central impact parameters in the range 0 < 6 < 0.75fm, 
with the purpose of restricting to the small portion of 
geometric cross section where the contribution of heavy 
residues is not dominant, and where IMF formation is en¬ 
hanced. The remaining fraction of cross section favour¬ 
ing compound-nucleus decays can be efficiently described 
through statistical approaches. Such a choice is how¬ 
ever schematic because, due to fluctuations in the cas¬ 
cade trajectories, the impact parameter is not directly 
characterising the entrance channel, and violent collisions 
may arise also in less central configurations with smaller 
probability. Conversely, less excited configurations are 
also associated to central impact parameters with smaller 
proportion than in peripheral collisions. A statistics of 
about 1500 stochastic events per system have been col¬ 
lected, using a 32 CPU parallel computing station. 



FIG. 5. Study of the fragment configuration for one event of 
the system ^®^Au-|-p at 750 A MeV, for reaction times which 
are compatible with the definition of the fragmentation pat¬ 
tern. The event shown in the figure is selected among those 
giving the largest fragment multiplicity. 


A. Dynamical description up to the formation of 
primary fragments 

Within the model described above, in fig. we repre¬ 
sent one possible evolution of the density profile of the 
systems ^®^Au-|-p at 750 A MeV for a central impact 
parameter; this is a rather rare event corresponding to 
the fragmentation of the target nucleus in more than 
three fragments. The system breaks into three asym¬ 
metric parts visible at 400 fm/c. At later times, further 
splits may proceed from some individual largely deformed 
sources, as displayed in fig. [^for the time 500 fm/c. In 
these spallation processes the fragment multiplicity sat¬ 
urates after 700 fm/c. This is shown in fig. (bottom), 
in correlation with the particle emission and the corre¬ 
sponding reduction of bound mass as a function of time 
for all the simulated systems (top). The middle panel of 
the figure shows the probability of observing a split in 
the system, as a function of time, for two of the reactions 
considered. 

The cooling process is reflected in the decrease of the 
average thermal excitation energy per nucleon shown in 
fig- HI A backbending appearing between around 50 and 



FIG. 6. Upper panel. As in fig. evolution of the mean 
number of emitted nucleons per interval of time for the spalla¬ 
tion reactions described in the text. Middle panel. Probabil¬ 
ity of split as a function of time for p-|-^°®Pb and p-|-^^®Xe at 
IGeV. Lower panel. Saturation of the fragment multiplicity 
Mfrag for p-|-^°®Pb at IGeV; Mean and standard deviation are 
shown for the multiplicity of fragment with Z>4. 
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FIG. 7. Evolution of the average excitation energy {E*/A) 
for the fraction of bound matter during the reaction. The 
double arrow gives the average uncertainty in terms of stan¬ 
dard deviation. The bunch of lines extending over the whole 
time evolution describes events where only a heavy residue is 
present; bunches of lines for systems which split into Muag 
equal to two, three or four fragments are divided by Muag 
for better visibility as they would all collapse on the line for 

Mfrag = l. 


100 fm/c indicates the attempt of the system to revert 
the initial pure expansion dynamics into the mechanism 
of fragment formation. Indeed, in presence of instabili¬ 
ties, it is energetically convenient for the system to break 
up into fragments. This also causes a slight increase of 
the temperature and thus of the thermal excitation en¬ 
ergy. Event by event, we consider as freeze-out time, 
tstop, the instant between t = 400 fm/c and t = 600 fm/c 
where the last split has occurred. Our choice is motivated 
by the fact that at t« 400 fm/c the split probability is 
maximum, whereas at twOOO fm/c it reduces to a quite 
low constant value. For events where a residue is ob¬ 
served, we adopt tstop=400 fm/c. Beyond the time tstop 
the decay process slows down and only sequential binary 
splits become possible, which can be efficiently described 
through a transition-state model. The dynamical calcu¬ 
lation is therefore completed with the model SIMON [74] , 
which incorporates in-flight Coulomb repulsion. 


B. On the way to the residue corridor 

In correlation with the excitation energy, also the 
isospin content of fragments and residues evolves in time. 
In general, when a compound nucleus is formed, its exci¬ 
tation energy is extinguished in an attempt of balancing 
proton and neutron decay widths, so that the bound mat¬ 
ter of the systems tends to accumulate along the residue 
corridor m, which is located in the neutron-deficient 
side of the nuclide chart with respect to beta stability, 
and any further decay occurs only along this line in av¬ 


erage. If however part of the excitation energy is spent 
in fragmenting the system, neutron rich fragments stop 
their decay path before reaching the residue corridor, in 
locations of the nuclide chart which are closer to beta 
stability, or that are even neutron rich m- 

This process inspired several experiments and simula¬ 
tions with statistical models where an assumption of ther¬ 
mal equilibrium of the system was imposed and a temper¬ 
ature was assigned [75] (the so-called ‘limiting tempera¬ 
ture for fragmentation’, corresponding to about 5 MeV). 
The dynamical approach handles this process without 
any hypothesis of equilibrium. Fig. (top left panel) ex¬ 
amines the evolution of the isotopic content for the six 
different hot systems for central impact parameters: the 
average isotopic content of bound matter, obtained by 
dividing the average number of bound neutrons consti¬ 
tuting the system fVbound by the average bound charge 
2^bound is tracked as a function of time until 700 fm/c. In 
this interval of time the path moves in average in the di¬ 
rection of the residue corridor while removing mass. The 
whole distribution of the isospin content {N) jZ of hot 
fragments is given in fig. for ^°®Pb-|-p, for the neutron- 
deficient system ^^^Xe-|-p and for the neutron-rich sys¬ 
tems ^^®Xe-|-p (top right, bottom right and bottom left 
panels); the following times are analysed: t = 200 fm/c, 
before fragmentation, t = 400 fm/c, after fragmentation, 
and t = 700 fm/c, when the fragment multiplicity sat¬ 
urates. In all the three systems, the distribution at 
400 fm/c covers the region of neutron-rich nuclei as a flat 
function of the element number and its distance from the 
residue corridor depends on the isospin content of the 
target nucleus; it drops to smaller values of {N)/Z for 
later times. As a function of the available excitation (i.e. 
of the time), the corresponding distribution of cold frag¬ 
ments is found displaced in the direction of the residue 
corridor. The cold IMF’s align along the residue cor¬ 
ridor only for the neutron poor target ^^^Xe; they do 
not reach it completely for the neutron rich target ^^®Xe, 
ending their decay path in the vicinity of /I-stability; in 
the case of the heavy neutron rich target ^°®Pb, only the 
largest fragments can reach the residue corridor at the 
end of their (shorter) decay path. A complete study of 
this reaction, involving also less violent events for periph¬ 
eral impact parameters would extend the distribution of 
residues to larger mass numbers which would accumulate 
along the residue corridor. The same behaviour charac¬ 
terises also the other heavy systems (not shown) and it 
recalls closely the experimental results for peripheral rel¬ 
ativistic heavy-ion collisions [ZillT]. 

We may also suggest that these results, in good 
agreement with previous studies based on statistical ap¬ 
proaches, indicate that the transport description is well 
adapted to follow the reaching of equilibrium conditions, 
through a chaotic population of the available phase space, 
within the dynamical process |78|. Up to this stage, this 
analysis agrees with inclusive data and statistical simu¬ 
lations, but it is not sufficient to characterise the mech¬ 
anism: both fission and multifragmentation can in fact 






FIG. 8. Top left. Evolution of the average isotopic con¬ 
tent of bound matter constituting different hot systems as 
a function of time (calculated for central impact parame¬ 
ters). Top right, bottom right, bottom left: Distributions 
of the average isotopic content of single elements produced in 
the systems ^^"‘Xe, ^®®Xe and ^°®Pb (moving clockwise) bom¬ 
barded by 1 GeV protons as a function of the element num¬ 
ber at 200fm/c (before fragmentation), at 400fm/c (latgest 
fragmentation probability), at 700fm/c (when the fragment 
multiplicity saturates in the dynamical calculations) and af¬ 
ter secondary-decay progressing from tstop- The /3-stability 
and the residue corridor are indicated (see text). Residues 
and IMF regions are also indicated. 


populate the neutron-rich side of the nuclide chart due 
to the curvature of the /3-stability valley. 


C. Fragmentation in few IMF 

From the analysis of the multiplicity of fragments with 
Z > A at 700 fm/c, studied for central impact parame¬ 
ters, we found that the lighter systems (Xe) prevalently 


recompact into one compound nucleus, or they undergo 
binary splits with about one order of magnitude smaller 
probability, and multiple splits are rare. The heavier 
systems, despite also displaying some tendency to re¬ 
compacting, are on the other hand dominated by binary 
splits, and ternary splits are also relevant. This analy¬ 
sis is presented in fig. The evolution of the fragment- 
multiplicity spectrum is also shown as a function of time: 
we observe that, even if density inhomogeneities arise at 
earlier times, the system starts separating into fragments 
rather late, at around 300 fm/c. 

An insight about the asymmetry of the splits is pro¬ 
posed in fig. by analysing the size correlation among 
the three heaviest IMF’s, of mass number Ai, A 2 and 
A 3 , produced in the same event, for events where at least 
two fragments are found in the range 5 < Z < 11. All 
combinations of the relative sizes fii, ^2 and ^ 3 , where 

= Ai/(Ai -I-A2 -I-A3), are used as coordinates in Dalitz 
plots. The size correlations are investigated for the sys¬ 
tems ^^®Xe and ^°®Pb, for central impact parameters, 
both for the hot (at 700fm/c) and for the cold systems. 
From this analysis we infer that, even when the frag¬ 
ment multiplicity is larger than two, the splits exhibit a 
large asymmetry. In the ^^®Xe hot system, represented 
by black contours positioned on the sides of the plot, at 
maximum two IMF’s are found in the range 5 < Z < 11, 




^frag (■Z>4) 


FIG. 9. Normalised yields as a function of the multiplicity of 
fragments with Z > A, for central impact parameters, at vari¬ 
ous time-steps of the dynamical process, and after secondary- 
decay. 
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5<Z<11 in >36Xe+p lAGeV 


5<Z<11 in 208pb+p lAGeV 


O 700fin/c 
• cold 


FIG. 10. Correlation among the three heaviest IMF’s of 
mass number Ai, A 2 and A 3 produced in the same event, for 
all events where at least two fragments are found in the range 
5 < ^ < 11. Combinations of relative sizes pi, p 2 and ps 
(where pi = Ai/[A-i-\-A 2 -\-A 3 )) are studied in a Dalitz plot for 
the systems ^^®Xe and ^°®Pb, for central impact parameters. 
Colour maps refer to cold systems after secondary-decay and 
the configurations at 700fm/c are indicated by black contour 
lines. 


and the fragment multiplicity is completed by a heavier 
residue. The action of the secondary decay may turn 
some few events into three-IMF patters which enter the 
selection and fill the centre of the plot. In the ^°®Pb 
system, hot and cold, symmetric splits are still rather 
rare with respect to events where one heavier fragment 
is present. The configuration of the splits has an obvious 
consequence on the kinematics. 


D. Charge distribution and kinematics: 
two emission modes for IMF’s 

The fragment-mass yields are shown in fig. The 
spectra at tstop and at the end of the sequential decay are 
similar except for the extremities, corresponding to the 
lightest and the heaviest masses, which have been mod¬ 
ified by a prominent light-particle evaporation process 
and by asymmetric fission. It is interesting to notice that 
the heavy-residue region is already filled at t = 400 fm/c, 
whereas IMF’s are also produced at later times. More¬ 
over, their final yield, after de-excitation has been con¬ 
sidered, is quite close to the yield given by the BLOB 
simulations at tstop- Therefore, within our calculation, 
the kinematics of the cold IMF’s should mostly reflect 
the kinematics of the hot IMF’s, when they are related to 
the most violent entrance channels. As shown in fig. [g 
the kinematics reveals therefore the explosive character 
of the process and is then modified by the Coulomb prop¬ 
agation. 

We conclude the analysis by recalling the initial inspir¬ 
ing experimental finding of fig. Due to the computa¬ 
tional complexity, we could not collect enough statistics 
to reproduce the same kinematic observable of fig. [l] for 
single isotopes, but we could produce a similar observable 
by collecting, for instance, the velocity distributions of all 


isotopes of carbon and fluorine for the system ^®®Xe-|-p at 
I A GeV, and all IMF’s with 7 < Z < II (interval chosen 
around oxygen and neon, which are elements frequently 
produced in multifragmentation) for the system ^®®Pb-|-p 
at I A GeV: this study is illustrated in fig. in the ref¬ 
erence of the heavy nucleus before the collision; the shift 
with respect to zero corresponds therefore to the mean 
recoil of the target. The spectra could be symmetrised 
because the global reaction configuration studied with a 
central impact parameter is symmetric. 

Garbon and fluorine in the ^®®Xe system change from 
a wide-hump distribution, for the hot IMF’s, to a con¬ 
cave distribution for the cold IMF’s. When present, the 
contribution of events with fragment multiplicity larger 
than two are indicated. For the ^®®Xe system it appears 
only in the cold system, leading to a wide convex por¬ 
tion of the spectrum, especially in the fluorine case: the 
convexity results from the variety of possible sizes and 
patterns involved in the splitting configurations, mainly 
when Mfrag > 2. The resulting overall concave or two¬ 
humped wide spectra of the cold carbon and fluorine iso¬ 
topes is produced by imparting different boosts to the 
fragments issued of binary events as a function of the 



FIG. 11. Production yields and kinetic energies for five spal¬ 
lation reactions, calculated as a function of the mass number, 
for central impact parameters, at 400fm/c, at the time tstop, 
and after secondary-decay (cold). 
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FIG. 12. Zero-angle invariant velocity distributions cri(v||)/CT 
of carbon and fluorine isotopes calculated for the system 
^®®Xe-l-p at 1 A GeV and of IMF’s with 7<Z<11 calculated 
for ^°®Pb-|-p at 1 A GeV. Left panels present distributions of 
hot fragments at 700 fm/c, while right panels present distri¬ 
butions of cold fragments after secondary decay. The integral 
of the distributions reflects the same number of events ex¬ 
ploited for the two systems and, for better comparison, they 
are scaled by the factors indicated in the boxes. The spectra 
are shifted with respect to zero by the mean recoil velocity 
of the ^^®Xe and ^°®Pb (indicated by arrows). The values of 
Mfrag indicate the multiplicity of fragments of Z > 4 (includ¬ 
ing heavy residues) associated to the events. The contribution 
for different fragment multiplicities is indicated. 


partner size, producing wide humps from the folding of 
different Coulomb boosts, and by an additional contribu¬ 
tion from asymmetric fission of the heavy residues, which 
selects a narrower Coulomb peak. From the analysis of 
fig. we infer that, even when the multiplicity is larger 
than two, the kinematics of the splits should however 
manifest a binary-like character due to the size asym¬ 
metry among fragments: the kinematics reflects in this 
case the prominent Coulomb repulsion imparted by the 
largest fragment. This effect becomes dominant in the 
208pb system, where concave wide spectra are also ob¬ 
served for larger fragment multiplicities. The calculation 
was limited to a small interval of impact parameters. The 
extension to the full range of impact parameters would, 


firstly, add or enhance the feeding of Coulomb peaks in 
the cold-IMF spectrum from asymmetric fission of heavy 
residues. Secondly, it would produce a folding over a span 
of recoil velocities for the target. Events where IMF are 
produced are related to a large range of central to semi¬ 
central impact parameters, and are mostly contributing 
to the centre of the distribution. Thus such folding would 
deform the central portion of the spectrum into an asym¬ 
metric shape with more extended tails for negative values 
of the velocity. In general, we observe that the more or 
less pronounced filling of the centre and the appearing of 
wide humps in the zero-angle spectra signs the presence 
of mechanisms possibly related to the sudden production 
of a few IMF’s in a same short interval of time, as sug¬ 
gested in experimental observations 


V. PHENOMENOLOGY 

Even though a quantitative comparison with experi¬ 
mental data is beyond the purpose of this work, we ob¬ 
serve that the model describes a large range of observable, 
from nuclide production to kinematics, which is globally 
consistent with the available experimental information. 

The dynamical evolution leads to a chaotic popula¬ 
tion of the available phase space, which makes the fi¬ 
nal result quite similar to the predictions of statistical 
multifragmentation models [7H] (statistical investigations 
along this line can be found in refs m)- The interest of 
the dynamical approach is that it can be used to extract 
further information on the phenomenology of the process 
at any time. Moreover, kinematical effects connected to 
the expansion dynamics can only be described within a 
dynamical model. 


A. Frustrated multifragmentation 


Fig. 13 gives an overview of the richness of the dy¬ 
namic behaviour. In the first instants which follow the 
collision, low-density tails appear in correspondence with 
the emission of preequilibrium particles, proceeding from 
forward angles and later extending more isotropically to 
all solid angles. 

After this time, the system starts expanding and the 
dynamical fluctuations handled by the BLOB treatment 
become a dominant mechanism in the process. With re¬ 
ducing bulk density, phase-space fluctuations grow in am¬ 
plitude and potential ripples develop, becoming the nest¬ 
ing sites of fragments: inhomogeneities in the bulk den¬ 
sity profile stand out at around 100 fm/c, but it takes 
them long time to eventually separate into fragments. 
The process exhibits a typical characteristic of the spin- 
odal instability, i.e. the arising of blobs of similar size in 
the bulk. The inset of fig. analyses for this same sys¬ 
tem the density averaged over those blobs as a function 
of their size at different times. The blobs are identified 
as any potential concavity found in the system and the 
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FIG. 13. Top. Time evolution of the density profile of 
the systems p (lGeV)+^°®Pb for one specific event selected 
among those giving multifragmentation. The system under¬ 
goes a spinodal behaviour, visible from 100 to 200 fm/c, when 
developing inhomogeneities of comparable size. Later on, the 
fragmentation mechanism is frustrated by the mean-field re¬ 
silience, resulting into a rather asymmetric fragmentation. 
Bottom. Time evolution of the size of potential concavities 
associated to the evolution of the density profile at the be¬ 
ginning of the process {t < 5 fm/c), during the phase of in¬ 
stability growth (t = 200 fm/c), and when fragments appear 
(f = 400 fm/c). See text for details. 


comparable size and preserve the spinodal signal also in 
the exit channel. We can observe that this is definitely 
not the case: not all blobs succeed in separating in single 
fragments but they bond together in groups. The event 
of fig, finally results in the fragmentation of the system 
into four asymmetric parts. Such a scenario seems to be 
general for this kind of spallation mechanisms and it is 
reflected in the low multiplicity of hot fragments analysed 
in f ig. and in the mass distribution of hot remnants of 
fig. |ll| In fact this latter, even when displaying a peak 
around the elements selected by the spinodal instability 
like oxygen and neon, presents a rather flat distribution 
which implies a large recombination of the spinodal in¬ 
homogeneities into larger fragments. As it was already 
argued in Ref. [59] . we can conclude on this phenomenol¬ 
ogy that the multifragmentation mechanism in spallation 
in the 1 AGeV energy regime is frustrated by the action 
of the mean field which tends to recompact the system as 
far as not enough energy is spent in the radial expansion. 


B. Exit-channel chaos and binary events 


Phase-space fluctuations favour fragmentation. More¬ 
over, they also act, over many events, in expanding the 
bundle of dynamical reaction paths into a large chaotic 
pattern of bifurcations [55]: this leads to a variety of 
exit channels. Fig. illustrates this effect by plotting 
the average radius of inhomogeneities found in the den¬ 
sity landscape as a function of time for several stochastic 
evolutions of the systems ^°®Pb-|-p at 1 A GeV. Each 
trajectory could lead to a different exit channel, char¬ 
acterised by a different kinematics and larger or smaller 
fragment multiplicities and charge asymmetries. 

Particularly interesting are trajectories leading to only 
two fragments in the exit channel, as it can occur rather 
often in systems like ^^®Xe-|-p at 1 A GeV, as illustrated 


size is their average radius (their shape is nearly spher¬ 
ical). At early times potential concavities coincide with 
the whole expanding system or with some large portions 
of it when particle flow develops. At late times poten¬ 
tial concavities have a large probability to coincide with 
the inhomogeneities arising in the density landscape: be¬ 
cause their size reflects the leading instability mode izni, 
they are all comparable in size, corresponding with larger 
probability to neon or oxygen nuclei [301157]. At inter¬ 
mediate times, sizes range from the whole system to the 
size of the spinodal undulations in the density landscape. 
This coexistence recalls phase-transition signals and cor¬ 
responding results at Fermi energies, where suitable ob¬ 
servables (such as the asymmetry between the charges 
of the two heaviest fragments produced in one collision 
event [50] or the size of the largest fragment produced in 
one event [81]) have been proposed. 

However, only in presence of a sufficiently large ra¬ 
dial expansion these blobs can separate into fragments of 



FIG. 14. A random selection of some reaction paths followed 
by the system ^°®Pb-|-p at 1 A GeV on the time-dependent po¬ 
tential landscape, represented by the average size of potential 
concavities as a function of time. 
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in fig. Also in this case, a spinodal process may acti¬ 
vate and immediately enter in competition with the ac¬ 
tion of the mean field which tends to reverse the frag¬ 
mentation pattern into a compact shape. Most of the 
times, this frustrating process results in one single com¬ 
pound nucleus. Very seldom the mean field succeeds only 
partially in coalescing the inhomogeneities of the density 
profile, and one or a group of those separates into a frag¬ 
ment and leaves the system. Depending on the stochastic 
configuration of the fragmenting system, the partial co¬ 
alescence could recompact the inhomogeneities in many 
combinations resulting in different asymmetries. 

To characterise the final exit channel, we have com¬ 
bined the dynamical process and the additional sec¬ 
ondary decay, which we simulated through a fission- 
evaporation afterburner. It turns out that two mecha¬ 
nisms may coexist in the spallation systems we exam¬ 
ined. One is a multifragmentation process frustrated by 
the coalescence effect of the mean field, which constrains 
the fragment multiplicity to small values, rarely equal to 
three or four, more often equal to two: this work focused 
mostly on this process. The other mechanism, related to 
the production of heavy residues, is also a binary split, 
but it coincides with asymmetric fission from a compound 
nucleus in the secondary decay process [3]- For different 
causes, both processes contribute to the same IMF pro¬ 
duction when regarding the charge distribution and even 
when regarding the isotopic content. However, these two 
contributions are fundamentally different because of two 
reasons. First: They are separate in time; multifrag¬ 
mentation develops within the short time of the dynam¬ 
ical collision process, while the second occurs during the 
longer time of the fission decay process. Second: from a 
dynamical point of view binary splits emerging from fis¬ 
sion or from multifragmentation are different: fission of a 
compound nucleus is a trajectory in a deformation land¬ 
scape which passes through the development of a neck; 
on the other hand, multifragmentation leads initially to 
a mottling topology in density space, possibly driven by 
spinodal instabilities, and then the expansion dynamics 
attempts to keep this topology, in competition with the 
antagonist tendency of reverting into a compact shape. 
This would lead to binary channels, similar to fission, but 
obtained through re-aggregation processes. The latter 
could generally have different kinematic features with re- 



FIG. 15. Time evolution of the density profile for the sys¬ 
tems '^^^Xe-l-p at 1 A GeV for one events resulting into an 
asymmetric binary split. 


spect to standard fission processes. Experimentally, this 
difference would be reflected in the velocity spectra of 
the IMF’s if the resolution is sufficiently high. 

The length of the transport calculation makes it pro¬ 
hibitive to track all contributions to the yields coming 
from the full distribution of excitation energies of the 
spallation system, and we had to restrict to the most vi¬ 
olent collisions only. More quantitative simulations are 
left for further works. 


VI. CONCLUDING REMARKS 

The Boltzmann-Langevin one body (BLOB) approach 
has been applied to the dynamics of hot nuclear systems 
produced in spallation reactions. A simplified cascade 
procedure was introduced to simulate the collision of rel¬ 
ativistic protons and light projectiles on heavy target nu¬ 
clei. Successively, the dynamics is followed in time within 
the BLOB treatment, to investigate the trajectory fol¬ 
lowed by the system as a function of the initial energy 
deposition. 

We find that unstable isoscalar modes can actually 
arise in such systems, like heavy nuclei bombarded by 
protons and deuterons in the I A GeV regime and, with 
a low but not negligible probability, become responsi¬ 
ble for the fragmentation of the system. According to 
our theoretical approach, we find that unstable modes 
in spallation should exhibit quite a similar phenomenol¬ 
ogy as spinodal instability in dissipative central heavy- 
ion collisions when the incident energy corresponds to the 
threshold between fusion and multifragmentation | 60 j . 

In the case of dissipative ion-ion collisions, the excita¬ 
tion of the system is mostly determined by mechanical 
perturbations, while in the spallation process the exci¬ 
tation originates from an almost isotropic propagation 
of the energy deposited by the light projectile. Though 
these situations are different in some aspects, they both 
drive phase-space fluctuations of large amplitude and 
they may both activate the spinodal behaviour and the 
amplification of mechanically unstable modes. In both 
cases, the mean field may then have the effect of revert¬ 
ing the whole system, or part of it, into a compact shape, 
smearing, modifying, or completely erasing the fragment 
configuration. In particular, if the kinetic energy feeding 
the expansion dynamics is not sufficient to disintegrate 
the system but it is still larger than the system can hold, 
the fragment multiplicity reduces and the fragment con¬ 
figuration becomes asymmetric. As an extreme situation, 
the process may look like asymmetric fission, but the 
chronology, as well as the violence of the process will be 
incompatible with the conventional fission picture: this 
process would correspond to a binary channel obtained 
by re-aggregation, keeping some dynamical aspects. 

Further binary contribution to IMF deriving from fis¬ 
sion have been accounted for by adding to the dynam¬ 
ical simulation a fission-evaporation treatment for the 
decay of all the possible warm configurations explored 
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by the collision. The overall simulation, constructed 
around a Boltzmann-Langevin description of the most 
violent collision events, leads to a correct qualitative 
picture for the production and kinematics of IMF nu¬ 
clides. In particular, the physical description that we sug¬ 
gest for IMF production in spallation may solve appar¬ 
ent discrepancies between some experimental interpreta¬ 
tions of inclusive and exclusive data, the former reveal¬ 
ing multifragmentation-like kinetic energies, the latter 
revealing small fragment multiplicities compatible with 
compound-nucleus decays. These experimental results 
would actually be perfectly coherent. 

As a conclusion, we found within a dynamical trans¬ 
port approach incorporating a Langevin term that the 


incident energy of a light projectile in the I GeV range is 
sufficient to turn the heavy target into an unstable sys¬ 
tem, where mechanically unstable modes develop lead¬ 
ing to a variety of aymmetric fragment configurations, 
including binary channels. This phenomenology reflects 
the entrance of the system in the spinodal zone of the nu¬ 
clear matter phase diagram and appears as a frustrated 
spinodal behaviour. 
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